Spontaneous spin polarization in doped semiconductor quantum wells 
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We calculate the critical density of the zero-temperature, first-order ferromagnetic phase transition 
in n-doped GaAs/AlGaAs quantum wells. We find that the existence of the ferromagnetic transition 
is dependent upon the choice of well width. We demonstrate rigorously that this dependence is 
governed by the interplay between different components of the exchange interaction and that there 
exists an upper limit for the well width beyond which there is no transition. We predict that some 
narrow quantum wells could exhibit this transition at electron densities lower than the ones that 
have been considered experimentally thus far. We use a screened Hartree-Fock approximation with 
a polarization-dependent effective mass, which is adjusted to match the critical density predicted 
by Monte Carlo calculations for the two-dimensional electron gas. 
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I. INTRODUCTION 

The interacting electron gas is one of the fundamen- 
tal systems of physics. However, in spite of a long tra- 
dition of study, the subject still has many open basic 
questions. Notably, the issue of the existence of a fer- 
romagnetic transition at low density has not been set- 
tled [3. Coulomb correlations play a central role in the 
low-density regime, and taking them into account the- 
oretically (i.e., going beyond Hartree-Fock) is unfortu- 
nately notoriously difficult. This problem has been most 
reliably tackled with numerically intensive Monte Carlo 
(MC) techniques H S 0. For the two-dimensional 
electron gas (2DEG), MC calculations indicate that, at 
T = 0, a first-order phase transition takes place at a cer- 
tain critical value r^^ of the dimensionless average sepa- 
ration between electrons = l/y/nNga'^, where Ng is 
the surface density and a|j is the effective Bohr radius in 
the embedding medium {a^ — 98. 7A for GaAs). 

The most widely used methods in MC calculations 6] 
are the variational Monte Carlo (VMC), which predicts 
a first-order phase transition at r^c — 13±2 {N^c — 
1.9 X lO^cm"^), and fixed-node diffusion Monte Carlo 
(FN-DMC) with which r^c = 25 (TV^c = 5.2 x 10** cm-^) 
has been found 0]. The VMC method uses a stochas- 
tic integration to evaluate the ground-state energy for 
a given trial wave function. The other method, which 
provides lower and more accurate ground-state energies, 
uses a projection technique to enhance the ground-state 
component of a trial wave function. In addition, in the 
FN-DMC method implemented in reference 0] , backflow 
correlations are included in the Slater determinant of 
the trial wave function, i.e. correlations are taken into 
account at the starting point of the process, for each po- 
larization. 

Of course, the ideal, purely two-dimensional electron 
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gas cannot exist in nature. Instead, experimentally 
one studies quasi-two-dimensional electron gases (quasi- 
2DEG) like the ones formed in modulation-doped semi- 
conductor quantum wells (QWs) Q. To the best of our 
knowledge, no Monte Carlo studies comparable to the 
ones mentioned above have been done for quasi-2DEG, 
but the ferromagnetic transition in QWs has been stud- 
ied theoretically in the frame of the local-spin-density 
approximation :9|. The critical densities predicted with 
that technique exceed by far the density interval given 
by MC for the 2DEG. 

On the experimental front, the spin susceptibility has 
recently been measured in GaAs/AlGaAs superlattices 
[Tol |. with electron densities as low as 1.7 x lO^cm^^. In 
spite of the fact that this value falls into the density range 
predicted for a transition by the 2DEG-MC calculations, 
no transition was observed. 

In this work, we study theoretically the possibility of 
a first-order transition at T = for the quasi-2DEG con- 
fined in GaAs/AlGaAs QWs as a function of the well 
width. We find that the width and the depth of the 
well play a crucial role in the existence of the transition. 
We prove rigorously the origin of the pronounced depen- 
dence of the transition density on the well width and 
predict the existence of an upper limit for this param- 
eter beyond which the polarized phase is energetically 
unfavorable. We find that the transition should happen 
at electron densities lower than those attained experi- 
mentally so far jinj |. for an optimum value of the well 
width, which we provide below. In our calculations we 
use a screened Hartree-Fock approximation scheme that 
includes a polarization-dependent effective mass which 
is introduced in order to take into account more accu- 
rately the effects of Coulomb correlation inside the well. 
This novel approach has the virtue of allowing us to use 
the results of the existing numerical Monte Carlo studies 
in two dimensions and extend them in a reliable way to 
quasi-2DEG systems. 

The paper is organized as follows. In Section FlI Al we 
introduce the basic scheme of the screened Hartree-Fock 
approximation and obtain the equations for the ground- 
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state energies for the 2DEG and the quasi-2DEG. In Sec- 
tion 111 Bl we discuss the results of this approximation. 
In Section IIIII we describe the polarization-dependent 
effective-mass approximation and we use it in combina- 
tion with the available 2DEG Monte Carlo data to make 
predictions for QWs. We end in Section Hvl with a sum- 
mary of our conclusions. 



II. SCREENED-HARTREE-FOCK THEORY 
WITH POLARIZATION-INDEPENDENT 
EFFECTIVE MASS 

A. Formalism 

In a quasi-2DEG, the HP equation may be written 

as [ni 
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where are the nth subband eigenstates and 

En \k) the corresponding eigenenergies, ml = O.OGTwe 
is the effective mass (me being the electron rest mass), 
e = 12.5 is the dielectric constant, e is the electron 
charge, A is the crystal area, and k is the in-plane wave 
vector. In all our calculations we take the z-axis as 
the growth direction of the heterostructure. The self- 
consistent potential vjP{z) is obtained by integration of 
the Poisson equation and is expressed as 
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where the ^-dependent electron density is 
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Here, A'';, is the doping sheet density and the external 
potential Vext{z) is the sum of the confinement poten- 
tial of the heterostructure plus the electrostatic poten- 
tial generated by the ionized donors (located symmetri- 
cally). The Permi level fc^^ for each subband satisfies 

2n{l + C)Ns = J2n ^Fn ■ The spin polarization index 
only takes the values C = and C = 1 since, at T = 0, no 
stable partially-polarized phases are possible in 2DEG 
0,^2- III contrast, recent calculations in 3DEG show 
that the transition is not of first order, but rather a con- 
tinuous one involving partial spin-polarization states @. 



Here we make the hypothesis that the quasi-2DEG be- 
haves like the 2DEG provided that the well width remains 
sufficiently small. 

We solve equation ^ following a method similar to 
that developed in reference llj. We expand the eigen- 



functions ^„i!{z) in the single-electron QW-basis func- 

(4) 



tions {0„(z), €„}, i.e. 

<i^w = E4"(^)^w- 



This enables us to write the eigenvalue equation 

H«'>{k)bi<Hk)^EiP{k)biP{k), (5) 

with b^Sik) = {a[^,hk), ■■■,a^pn{k), ...f- The matrix ele- 
ments of the Hamiltonian operator H^'^\k) are 
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We note that equations from reference lll| have a num- 
ber of misprints which are corrected here [l3|. 

To introduce screening in the HP approximation we 
note from equations Q) and © that Vfp^^'^\k) is the 
only term present in the pure 2D case and describes the 
long-range in-plane Coulomb interaction. Thus, we have 
dressed the interaction line of the exchange diagram 
by replacing the bare Coulomb potential V{q = |k — 
k'l) = {2TTe'^ /e){l/q) in equation JSJl with the statically 
screened Coulomb potential 
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where q. 



(C) 



(2 — C,)/a*g is the ^-dependent Thomas- 



Permi wave number for the 2DEG. Unfortunately, the 
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FIG. 1: Different Coulomb potentials. The solid line corre- 
sponds to the unscreened Coulomb potential V{r) = 1/r. The 
dashed (dot-dashed) curve corresponds to the Thomas-Fermi 
screened Coulomb potential for the polarized (unpolarized) 
case. The dotted curve represents the potential V{r,dw), 
Fourier transform of equation (with respect to g = |k — k'|) 
for an infinite QW of dw = lOOA, for t = r = q = p= 1. 
Inset: the Fourier transform of equation @ for an infinite 
QW of dw ~ lOOA for different values of t, r, q,p. 



Thomas-Fermi screened Coulomb potential for the po- 
larized (dashed hne) and unpolarized (dot-dashed line) 
cases, for t — r = q= p= l. In the inset we show that 
this case dominates over all others (we only show a few 
relevant examples). 

In what follows we assume that only the first subband 
is occupied, and therefore the summations over the sub- 
band index may be omitted. 

In solving the eigenvalue equation by iteration, we con- 
sider that self-consistency is achieved at the Zth step when 
\ai^!^'\k) - 4i^^'"'^(fc)|/|4i^^'-')(fc)| < 10-4 for all p,n 



and k. We obtain the set of eigenfunctions and eigenener- 
gies {^l^il{z),EP{k)} for C = and C = 1, which allows 
us to write the ground-state energy per particle 
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Fourier transform in 2D real space of equation H1U|) can- 
not be obtained analytically but its decay at large r is 
found to be [l3 
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On the other hand, Vt^p^'-'^\k) (Eq. 0) which arises 
from the intrinsic inhomogeneity of charge distribution in 
a quasi-2DEG, represents an interaction of short range. 
To see this, we consider the Fourier transform in 2D real 
space 
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of the potential V{q,z- z') = (27reV£)(e-9l^-^'l - \)/q 
contained in equation (|5Jl. It can be shown that V{r^\z — 
z'\) is of short range in the plane |3| and that at large r 
it can be approximated as 



=2 j2 



V{r,dw) = -0.032^^ — 



(13) 



where we calculate the z and z' integration of \z — z'\'^ 
for an infinite QW of well width dw setting the sub- 
band indexes equal to one. From equations (|ll|l and 
(|13|l we note that both potentials decay at the same 

rate at large r and that V{r,dw) ~ 0.128Vs^\r) and 

V{r,dw) ~ 0m2V}^\r) for dw = a*^. 

We show in Figure that the potential defined by 
equation ((T^ . suitably integrated over z and z' (dotted 
line), has a range that is shorter than the range of the 



In order to check our quasi-2DEG calculations in QWs 
when dw tends to zero, we shall compare our results to 
the 2DEG screened HF case [13 
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1±^-— [2C + V2(l-C)]I(:rc) 
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where a:,^ = ilC + 2\/2(l — C)]^s is the polarization- 
dependent Thomas-Fermi wave number divided by 2fc^^ 
and 

f ^ xdx I 

I(^c) = / [arccos(x) — 1 — a:;^]. (17) 



B. Results 

In this section we present the results obtained with the 
screened Hartree-Fock approximation with polarization- 
independent effective masses introduced in the previous 
section. 

In Figure|21 we plot the critical density for infinite QWs 
(in which there is no barrier penetration) in the screened 
HF approximation as a function of well width dw (solid 
squares). The limiting point at dw = (2DEG) was 
calculated with equation (|16|l and the remaining points 
with equation (|14() . The excellent match between these 
two different equations reflects the correctness of our 
derivations and calculations 13]. We see that the criti- 
cal density for the infinite QWs is a rapidly decreasing 
function of dw ■ This result, which looks superficially 
correct on the basis of a gradual 2D to 3D transition, 
can be explained rigorously in terms of the interplay be- 
tween different components of the exchange interaction 
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FIG. 2: Well- width dependence of the critical density in 
the screened HF approximation. Solid squares correspond 
to infinite quantum wells and solid circles to finite QWs with 
well height 14=247 meV. Open squares correspond to infinite 
QWs with polarization-dependent effective masses with a con- 
stant ratio / = ml/rriQ = 0.65. The solid and dashed lines 
are obtained with equations 1241 and 12711 (with / = 0.65), 
respectively, for infinite QWs. 



in a quasi-2DEG, as follows. We found that the coef- 
ficients a^j^-^{k) (Eq. |0J) change very little with k and 



a^i{k) « 1 whereas a'^i{k) <C 1 for p > 1. This result is 
a direct consequence of the small-density condition since, 
if this condition is met, the ground-state wave function 
must retain the shape of its one-electron counterpart in 
a QW. Thus, V}^''^'^\k) (Eq. ®) is positive and it is the 
leading matrix element, yielding a negative (see Eq. ©) 
and dvi'-independent contribution to the exchange en- 
ergy (pure 2D case). In contrast, the matrix elements 

v}p'^^''\k) are slowly varying functions of k since we are 
studying small-width and low-density heterostructures. 
If these conditions are fulfilled, we can expand the expo- 
nential in equation up to first order since the exponent 
satisfies the following inequality: 



|k - k'||z - z'\ < 2k^pdw = 2^2(1 + 



< 0.3, 



(18) 

provided that > 13 (lower MC limit for the transition) 
and the well width dw < {— 98. TA). Then we per- 
form the integration over z and z' {infinite QWs) yielding 

Gil, 11 — — 0.207diy. Thus, V-j^i^^^^ is negative, propor- 
tional to dw and fc-independent. The matrix elements 



(i)(C) 



tp 



for indices with opposite parity are always zero 
since due to the symmetry of the wave functions we have 
Gti.ip — 0. On the other hand, we have evaluated Gti^ip 
for indices with equal parity obtaining that these are con- 
siderably smaller than Gn^n rendering V^^^^'-^'' practi- 
cally diagonal. Also, {(j)t\vjP \<pp) « (see Eq. ®) for 
the non-diagonal elements since, due to the low-density 
condition, v}p must be a very slowly varying function 
of z. Thus, H[p\k) is almost diagonal, being its first 
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+ (^i|yif)|0i)-^«(^)-yA'»'^)(fc). 

(19) 

With this expression for E^^\k) we can obtain an ap- 
proximate equation for Ngc- We only need to make two 
easily justified additional approximations. From the be- 
havior of the coefficients a:^^{k), i.e. a^^^''(fc) « aii''(fc) ~ 



1 and Opiik) 



(0) 



(fc) « for p > 1 it can be seen 
(from Eq. (HHJl) that e^P {k) « ei and (from Eq. ^) 
^^P{z) « 4>i{z)- Using the latter in equation Q we 
get n«)(2) « (2 - 0(1 + C)iVs</'?(^)/2 = N,(t>l{z), and 
therefore a ^independent {(f)i\Vsc\(j)i) (see Eq. Q). By 
inserting equation H19|l in equation (|14|l we obtain, after 
some algebra, the following equation for the energy shift 
between both phases 
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where 
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with 
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(22) 
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(23) 



where x = k/kp"^ in 50 and x = k/kp' in gi. The same 
holds for x' and k'. 

Taking into account that the energy shift E^jpp — E^^pp 
must be zero at the transition density, we may write the 
following equation that relates dw and Nsc 



1 + 0.207— = 4:F{N„ 



(24) 



Finally, to demonstrate that the transition density is a 
decreasing function of the well width, we need to prove 
that F{Nsc) is also a decreasing function. To see this, 
we note firstly that gi > g^ for all values of its argu- 
ments making the function F{Nsc) always positive allow- 
ing equation l|24|) to be solvable. Secondly, both go and 
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gi are decreasing functions of Nsc for all values of x', 
and it is straightforward to prove that gi decreases more 
rapidly than g^ making F{Nsc) a decreasing function. 
Thus, an increase of dw must be accompanied with a de- 
crease of Nsc proving that the monotonically decreasing 
dependence of Nsc on dw is governed by the competing 
action of the different components of the exchange inter- 
action: the in-plane component represented by F{Nsc) 
and the out-of-plane term driven by dw- The behavior 
of go, gi and F that we described can also be recog- 
nized in the unscreened HF case indicating that the de- 
pendence of Nsc on dw is purely due to exchange. The 
validity of our approximation can be verified in Figure 
El where the solutions of equation (|24|l are depicted with 
the solid line. This curve approaches very well the exact 
values (solid squares) indicating that the approximations 
we made to derive it are well justified. For low densi- 
ties (high Ts) and low well widths, the solid curve fits 
excellently the solid squares since in these regimes the 
approximations we made become exact. In the interme- 
diate region {dw ~ lOOA) the solid curve fits very well 
the solid squares. 

A look at equation H24|) also indicates that there exists 
an upper limit for the well width. In fact, due to the 
decrease of F{Nsc) and since F{0) is finite, it can be seen 
from equations 122|l and l|23|) that the following relation 
holds 



0.207- 



IWL 



(0) 

S 
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Thus, it must be dw < dwL — 2.42a^ « 239A to allow 
the confined electrons to reach the polarized phase. This 
new result is entirely due to the Thomas-Fermi screening 
and is not present in the unscreened HF approximation 
since F{0) diverges in this case. Furthermore, the ratio of 
the Thomas-Fermi wave numbers for both polarizations 
must be qf ^/qi^^ > 3/2 (in our case it is qi'^'^ /qi^^ = 
2). Otherwise, no polarized state could be possible in 
a QW. These key results indicate that the well width 
plays a crucial role in the search for spontaneous spin 
polarization in QWs [l8l |. 

The implications of the existence, according to our cal- 
culations, of an upper limit for the well width must be 
considered with some care. If one attempts to reach a 
3DEG system by increasing the well width one would ap- 
parently fall into the paradox that no transition is possi- 
ble in 3DEG. This conclusion is incorrect for two reasons. 
First, we must take into account that the 3DEG-MC re- 
sults are obtained in the jellium model, in which the pos- 
itive background is taken to be a uniform neutralizing 
static charge distribution, whereas in our quantum-well 
calculations the positive charges of the ionized donors 
are located far away from the electron gas, which re- 
sults in an important change in the direct Coulomb en- 
ergy. In other words, wide-enough quantum wells and 
3DEG jellium model must be considered as different sys- 
tems. Secondly, our calculation assumes that only one 



subband is occupied (a valid assumption in narrow quan- 
tum wells at low density) whereas any extrapolation of 
our conclusions to 3DEG systems would have to contem- 
plate necessarily occupation of many subbands. Thus, 
we reach the conclusion that the most likely scenario is 
that there is a phase transition in narrow quantum wells, 
which disappears for intermediate well widths, and reen- 
ters at wider well widths as expected when the 3DEG 
limit is approached. 

We now use the previous result to analyze the criti- 
cal density for finite QWs, plotted in Figure|21with solid 
circles. We set the height of the QWs to Vt = 247 meV, 
a typical experimental value [Toj. This curve exhibits 
a non-monotonic dependence on the well-width show- 
ing a maximum for dw ~ 35A. Also we observe a gen- 
eral reduction of the critical density with respect to the 
case of infinite QWs. This can be simply understood 
in terms of the previous result (monotonically decreas- 
ing critical density for the infinite wells) and the pen- 
etration of the electron wave function into the AlGaAs 
barriers; the latter causes the wave function to spread 
beyond the nominal well width, effectively "enlarging" 
the well. As a consequence, for example, a finite QW 
of dw ~ 60A has the same critical density as that of 
an infinite QW of dw ~ lOOA. In fact, the penetra- 
tion depth ds = hU/2mt {Vh - E^ ) increases when Ei 
is raised as dw is lowered [l6|. This effect produces an 
inflection point at dw ~ 75A and the mentioned maxi- 
mum at dw ~ 35A due to the competition between dw 
and ds- 

Let us go back to the curve for infinite QWs in Figure 
El The limiting {dw = 0) value Nsc = 17.5 x lO^cm'^ 
corresponds to Vsc — 4.32, showing a sizable increase with 
respect to the (unscreened) HF value Vsc — 2.01 [20l| . 
This increase, however, is not sufficient if we consider the 
value Tsc = 13 obtained in reference using VMC. This 
indicates that a significant degree of Coulomb correlation 
is being left out in the screened HF approximation. 



III. POLARIZATION-DEPENDENT 
EFFECTIVE MASSES 

A. Two-dimensional case 

In order to go further and improve our treatment 
of Coulomb correlation, we need an approximation 
scheme applicable to the quasi-2DEG such that as dw 
tends to zero (pure 2DEG) the critical density ap- 
Toaches the values predicted by MC calculations 0, 
~|]. To achieve this, we incorporate phenomenolog- 
ical polarization-dependent effective masses TOq (unpo- 
larized) and (polarized) in our formalism. Due to 
the lack of experimental data on effective masses in 
GaAs/ AlGaAs heterostructures for both polarizations 
and that no calculations on polarized effective masses 
exist in 2DEG, we resort to calculations of unpolarized 
effective masses and ground-state energies in pure 2DEG 
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FIG. 3: Ground-state energy shift Eshift ~ E^hf2d " Ehf2d 
versus Ts- Solid circles represent the reported values in Table 
I of reference |3. The error bars denote the VMC standard 
errors. Open squares correspond to the values tabulated in 
Tables I and II for the VMC method in reference Q . No error 
bars are plotted for clarity. Open circles belong from reference 
0. The solid and dot-dashed lines represent our calculations 
obtained with equation 12611 for polarization-dependent effec- 
tive masses using / = 0.65 and / = 0.49 respectively. The 
dashed curve corresponds to the same calculations but us- 
ing the values of / that come from /(r^) showed in the inset 
(dashed line). In the inset, solid and dashed lines correspond 
to the values of / that fit the curves Eshi ft from reference 
and reference 3 respectively. 

13, 0, 112 to justify this procedure. Let us summarize 
the conclusions of those studies relevant in our context: 

(a) Coulomb correlation increases the effective mass [2^. 

(b) The absolute value of the correlation energy of the 
unpolarized 2DEG ground state is greater than its polar- 
ized counterpart 0- 

(c) The absolute value of the correlation energy is greater 
in 2D than in 3D (both unpolarized), leading to 2D ef- 
fective masses substantially larger than those of the 3D 
case at equal plj . 

(d) The correlation-energy shift between both phases in 
2D is greater than the unpolarized correlation energy 
shift between 2D and 3D [||. 

Making use of (a) and (b), with the supporting evidence 
of (c) and (d) , we conclude that TOq must he substantially 
larger than m\ at equal . 

By defining the ratio / = m\/m'^ and rewriting xq = 
\[fC, + 2\/2(l — C)]rs we may write equation IjlBfl as 

(0 e^mg f 1 + C 

''^^^^-2a|,e\[(/-l)C+lH 

- — [2C + %/2(l-C)]I(a:c)V (26) 

We observe from reference 21;] that TOq ~ 1.2 for 
> 5 in the modified Hubbard approximation. That 
approximation is an attempt at including correlation ef- 
fects by means of the introduction of the Thomas-Fermi 
wave number in the so-called local-field correction factor. 



Since we have incorporated screening correlations and HF 
effects within a similar scheme, we take mg = 1 in equa- 
tion l|26|l to avoid an overestimation of the effective mass 
in the unpolarized phase. 

Using equation the lowest 2DEG-VMC value, i.e. 
rsc = 13, is obtained with / = 0.65. With this value of / 

we calculate Esm ft — £'lfF2D " -^hf2£) ^'^'^ P^*-"* versus 
Ts in Figure|21 (solid line). Here we have taken e^ /2a*^e = 
1 Ry to compare our Eg^i ft against MC results. The solid 
circles correspond to the results obtained in reference 
and the open squares are from reference Q. We note 
the excellent agreement between our curve and the MC 
points: by adjusting only one point our curve meets all 
the points obtained in reference 0, 0] ■ This agreement 
implies that the ratio between both effective masses de- 
pends weakly on the density and supports the validity of 
our assumption of a density-independent / factor. This 
conclusion is consistent with the fact that in the modified 
Hubbard approximation the unpolarized effective mass is 
a slowly varying function for Ts > 5 [IJ . If this were also 
the behavior of the polarized effective mass, we could 
conclude that / would be a slowly varying function of Ts . 

We now repeat the previous analysis but using the data 
of Attaccalite et al. [j| . In that paper the authors obtain 
rsc = 25 which, as we mentioned in the introduction, is 
the highest value found in the literature for spontaneous 
spin polarization in 2DEG at zero temperature. We find 
that equation (|26|l reproduces the value rsc = 25 when 
/ = 0.49. In Figure O (dash-dotted line) we plot the en- 
ergy shift, Eshift = e'^hf^d - Ehf2D^ versus r^, for this 
value of /. This curve does not fit well the data from 
reference shown as open circles. Instead, we find that 
the ratio / now taken as a function of (dashed line in 
the inset) fits very well the MC points calculated in ref- 
erence 3 (open circles) when it is used in equation H26() 
(dashed line in Fig. O For completeness, we show in the 
inset (solid line) the values of / as a function of rg that fit 
the parametrization of Eshift obtained by Ceperley 
This curve exhibits a weak dependence on rs for > 10 
giving support to our initial assumption of a constant 
value / — 0.65. 



B. Quasi- two-dimensional case 

We now apply the polarization-dependent effective- 
mass approximation to the quasi-2DEG by means of a 
slight modification in equation ^ . We note that the two 
effective masses ml on both sides of equation Q belong 
to different situations |[ij|: the ml on the l.h.s. represents 
the in-plane effective mass and therefore is being affected 
by Coulomb correlations. In contrast, the on the 
r.h.s. reflects the out-of-plane effective mass of one elec- 
tron moving in the z-direction governed mainly by Vext (z) 
and v}p{z), and thus not being affected by Coulomb 
correlations, according to our discussion about screening 
m Section lira Then we solve the eigenvalue equation, 
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FIG. 4: Well-width dependence of the critical density in the 
screened HF approximation with polarization-dependent ef- 
fective masses. Open squares (circles) correspond to infinite 
(finite) QWs with a constant ratio / = ml/mQ — 0.65. Up 
(down) triangles correspond to infinite (finite) QWs for a ra- 
tio / which depends on Vs (dashed line in the inset of Fig. 



equation and use equation (|14|1 for both polariza- 
tions incorporating the effective masses ttlq — ml and 
TO* = fml (the last one only in the in-plane terms). We 
plot in Figure [51 with open squares, the results obtained 
for infinite QWs for / — 0.65. We calculate the limiting 
point at dw — with equation 126|l and the remaining 
points with equation l|14|) with the above-mentioned re- 
placement giving an excellent match. We show in the 
Figure 01 in a logarithmic scale for the vertical axis, the 
results for finite QWs (open circles) and infinite QWs 
(open squares) where we have taken / = 0.65. Both 
curves exhibit the same general characteristics as in Fig- 
ure [3 (sohd squares and solid circles). 

Now we obtain an equivalent of equation (|24ll by in- 
corporating the ratio / that multiplies the band mass 
for the polarized phase, yielding 

2 1 + 0.207^ =4F(7V,„/), (27) 
/ 

where F{Nsc, f) is the same as in equation ((21]), but gi 
now reads 



gi(Nsc,x,x',(p) = \/TrN, 



iscfJ-B 



X y/(x + x')^ - Axx'sin'^ip + Jq^P /qf^ 



(28) 



The polarization-dependent effective-mass approxima- 
tion does not change the previous result regarding the 
monotonically decreasing dependence of Nsc on dw since 
/ < 1. We show that the solutions of the approximate 
equation (|77|l (dashed curve) depicted in Figure [21 for 
/ — 0.65, fit perfectly the exact values (open squares). 
We observe that since F(0, /) — j(j — |), dwL does 
not depend on / (see Eq. (|77|) ). Thus, dwL depends on 



correlations, in our Thomas-Fermi model, via the ratio 
qi^^ /qi^^ and, consequently, the relation dwL = 2.42a|j is 
a universal one, i.e. it holds for QWs of any material. We 
note that this interesting result does not depend on the 
approximations we made to derive equation (|27() since 
those become exact as Ns tends to zero. 

On the other hand, the location of the observed maxi- 
mum remains unchanged with respect to the / = 1 case 
(solid circles in Fig. ^ in accordance to our initial as- 
sumption that screening and correlations manifest only 
in the plane. Up (down) triangles correspond to infi- 
nite (finite) QWs where we have taken the ratio / as 
the dashed curve in the inset of Figure [3 (Attaccalite et 
al. 13]). We observe a drastic diminution of the transition 
densities in this case for both infinite and finite QWs. We 
note that the MC density interval for spontaneous spin 
polarization mentioned in Section ^ appears notoriously 
shrunk for the finite QWs studied here. In fact, from Fig- 
ure^ we obtain a new density interval for the transition 
densities in finite QWs between Nsc — 3.2 x 10''cm~^ 
and Nsc = 6.1 x lO^cm"^. We take these values from 
the transition densities at the maximum of the curves 
related to FN-DMC (down triangles) and VMC (open 
circles) respectively. 

In reference |10| . the spin susceptibility (= m*g*) has 
been measured in a high quality 200-fold GaAs/AlGaAs 
superlattice of 100 A of GaAs wells and 30 A barriers of 
Alo.32Gao.68As, with unprecedented low densities such as 
Ns = 1.7xl0^cm^^ (vs = 13.9) and no transition was ob- 
served. According to what we have mentioned above, this 
is not surprising. There are several possible reasons for 
this negative result. We first note that the density used, 
although low enough for a transition in the pure 2DEG, 
is clearly too high considering the finite well width for 
finite QWs: for dw = lOOA (Fig. [J), the elect ron den- 
sity achieved in reference (ICl] is 6.5 times higher than our 
critical density (open circles) which uses the ratio / that 
matches the 2DEG value from reference and 50 times 
higher than the critical density (down triangles) which 
uses the r^-dependent ratio / that matches the 2DEG en- 
ergy shifts Eshift from reference Also, due to the tun- 
neling of the electrons into the AlGaAs barriers, the su- 
perlattice acts like a single, extremely wide QW. Further- 
more, it is possible that if the QW were sufficiently wide, 
the quasi-2DEG could lose its two-dimensional character- 
istics, allowing for stable partially-polarized phases like 
those possible in the 3DEG, turning more difficult the 
detection of the transition. We note that the effects of 
in-plane correlations combined with the finite well widths 
and heights of QWs produce a drastic diminution of the 
transition densities by a factor that ranges from 3 to 15 
depending on which method VMC or FN-DMC turns out 
to be the best tool to estimate the transition density in 
pure 2DEG. For the best case, it should become nec- 
essary to achieve electron densities lower that the ones 
studied experimentally thus far by a factor of 3 and by 
a factor of 53 in the worst case. Very different could 
be the quasi-two-dimensional hole gas (quasi-2DIIG) see- 
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nario since in that system, high values such as « 80 
are aheady attainable . However, our theoretical pre- 
dictions about the critical transition density in n-doped 
GaAs/AlGaAs QWs cannot be straightforwardly trans- 
lated to quasi-2DHG. The adaptation of our formalism 
to the problem with holes is currently in progress. 

Based on the insight gained from our calculations, we 
propose that the optimal conditions for observing a fer- 
romagnetic transition in multiple QWs are: 

(a) well widths between 30A and 50A 

(b) wide AlGaAs barriers between wells to prevent tun- 
neling, and 

(c) well height VJ, as large as possible to minimize barrier- 
penetration effects. 

In a recent experimental work. Gosh et al. [23 | 
report a possible spontaneous spin polarization in 
mesoscopic two-dimensional systems that is at odds 
with our findings. They have used 2DEGs in Si (5-doped 
GaAs/ AlGaAs heterostructures with densities as low 
as Ns — 5 X lO^cm^^ {rg — 7.6) and the temperature 
was set at T 40 mK or equivalently T/Tp « 0.02 
since Tp ~ 2.3 K at rg = 7.6. Somewhat surprisingly, 
according to their interpretation of the data, these 
authors found partial spin polarization with C = 0.2. 
The authors attribute this partial spin polarization to 
the finite T since no partial spin polarization is possible 
in 2DEG at T = [1 However, in reference 

the authors find partial spin polarization for T/Tp 
between 0.3 and 1.6, i.e. well above T/Tp — 0.02 
reported in reference |24| . On the other hand, = 7.6 
is considerably lower than the lowest value for spin 
polarization in 2DEG . 



IV. SUMMARY 

In summary, we have calculated the ferromagnetic 
critical density at T = of the quasi-two-dimensional 



electron gas confined in semiconductor GaAs-based 
symmetrically-doped quantum wells. We use the 
screened Hartree-Fock approximation and prove rigor- 
ously that the pronounced decrease of the transition den- 
sity with the well width is governed by the interplay be- 
tween the in-plane and the out-of-plane components of 
the exchange interaction. The combination of these ex- 
change terms with the Thomas- Fermi screening produces 
a universal upper value for the well width beyond which 
the polarized state cannot exist. We add different effec- 
tive masses for both spin polarizations, which are intro- 
duced in order to take into account Coulomb correlations 
beyond screening. Once the value of the effective mass 
for the polarized phase is adjusted so as to reproduce 
the transition density for the pure 2D case calculated 
with the VMC method, our theory gives ground-state 
energy shifts that agree with those calculated within this 
method. On the other hand, a density-dependent ra- 
tio between both effective masses is required to fit the 
ground-state energy shifts calculated with the FN-DMC 
method. Based on our theory and the existing MC cal- 
culations for the 2DEG, we predict that narrow quantum 
wells (with well widths roughly in the range 30A< dw ^ 
50A) should exhibit a ferromagnetic transition at a den- 
sity range between Nsc — 3.2 x lO^cm^^ (r^ w 100) and 
Ngc = 6.1 X lO^cm"^ {ts ~ 23). This range, which looks 
far from the densities achievable nowadays in GaAs quasi- 
2DEG, is already within reach in GaAs quasi-2DHG sys- 
tems. 
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